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We investigate groundstate energies and magnetization curves in the one dimensional XXZ-model 
with next to nearest neighbour coupling a > and anisotropy A( — l<A<l)atr = 0. In between 
the familiar ferro- and antiferromagnetic phase we find a transition region - called metamagnetic 
phase - where the magnetization curve is discontinuous at a critical field Bc{a, A). 

on" 

a^ 

Q^ , I. INTRODUCTION 

T-H . 

,^ . Experimental results for the magnetization curves of FcxMni-xTiOs, GdNi2Sb2, GdCu2Sb2 or Tbi^xScxMn2 

1^ show a rapid increase (or discontinuity) if the applied S-field exceeds a critical value Be- IfQij-fB > B^ the substance is 

1"^ , almost fully magnetized. This phenomenon is called 'spin-flip' or 'metamagnetic' transitionJu. There have been made 

(^ ' various attempts to explain the 'metamagnetic' transition in the context of Ising-like Hamiltonians. It is the purpose 

of this paper to show that discontinuities in the magnetization curve can be seen as well in the one dimensional spin- 2 

XXZ-model with next to nearest neighbour (nnn) coupling 



> 

in 
in 

o 



N 
i=l 

O ■ ^h- {SfSf^2 + SfSl, + A5f 5f+2) + BS^ (1.1) 

0^. 

^\ . in the presence of a uniform external field B. We chose the next nearest (nn) neighbour coupling Ji to be antiferro- 
magnetic (Ji > 0) and use the notation a = J2/ Ji- In the a — A - plane, wc will primarily concentrate on the regime 

?5 : a> 0, -1 < A< 1. pp 

The isotropic model with A = 1 and nnn coupling a has been investigated by many authorsoS. Most of these 
investigations focussed on the transitional from the 'spinfluid phase' a <,«c-ito the dimerphase a > ac- The transition 
point Uc = 0.2411... has recently been determined with high precisiorJallJ by means of conformal field theory and 

Q ' renormalization group techniques. . .■ . 

Q , The Hamiltonian with a = 1/2 , A = 1 has been studied first by Majumdar and Goshll2Hl3. They found that the 



X3 



X 



'dimer states' 



\i^>=^[hmA]-[N~l,N] (1.2) 

^ |0>=^[2,3][4,5]...[iV,l] (1.3) 

are eigenstates of H{a = 1/2, A = 1, B ~ 0). Here 

[^,i + 1] = ^[x+Wx-(* + 1) - X-«X+(« + 1)] (1.4) 

are nearest neighbour (nn) valence bond states with total spin zero, called dimers. Van den Broektll proved that the 
dimer states are indeed g roun dstates of the Hamiltonian at the 'Majumdar-Gosh' point (a = 1/2, A = 1). Affleck, 
Kennedy, Lieb and TasakiE3ll3 were able to show that the dimerstates are the only ground states and that there is a 
finite gap to the first excited state. |— ■ 

Hamada, Kane, Nakagawa and Natsumellll discussed uniformly distributed resonating valence bonds (UDRVB) in 
the generalized railroad trestle model, which is equivalent to the isotropic linear Heisenberg chain with nn and nnn 
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interactions. They found that for negative Ji and J2 — —1/AJi the UDRVB is the ground state which is degenerate 
with the fully magnetized state with total spin S = Sz = N/2. As we will show later this phe nomenon also occurs 



for positive values of Ji if the parameters a and A are properly chosen in the Hamiltonian (1.1) 



Shastry and Sutherland discussed the frustrated model with differen interaction strengths in x,y and z directionll3. 
The critical propcrticsxtf the anisotropic model (A 7^ 1) in the absence of an external field B have been elaborated by 
Nomura and OkamotoE3. They confirmed-|that this model and the qu antu m sine-Gordon model belong to the same 



universality class. Tonegawa and HaradacJ have studied Hamiltonian (1.1) with ferromagnetic nn and antifcrromag- 
netic nnn interactions for positive A. 



The dimerstates ( |1.2D , (1.3) are eigenstates of the anisotropic model along the whole line a=i,— oo<A<ooas 
will be shown explicitly in section 2. 
However, the eigenvalues 

Eo{a^^,A)^-Ni^ + j) (1.5) 

are groundstate energies only for A > — i. For A < — i the groundstate |F± > with energy 

i?f(a = i,A) = -^AiV (1.6) 



is found in the ferromagnetic sector where all spins are down or up. This is a first hint, that the model (1.1) is 
particularly suited to study the transition from antiferromagnetism to ferromagnetism. 

The outline of the paper is as follows: In section 2 we report on the quantum numbers and the finite size effects of 
the groundstates as they depend on a, A a nd the magnetization M = S^/N . Section 3 is devoted to an analysis of 



the magnetic properties of the model (IT). Three phases can be found in the a — A-plane: the ferromagnetic, the 
antiferromagnetic and the metamagnetic phase. 

The Hamiltonians with a = 0, A = — 1 and a = 0.5, A = —0.5 are special in the sense that the groundstate is highly 
degenerate — namely with respect to Sz = 0,±1,±2, ...,±N/2. This feature is discussed in section 4. 

II. QUANTUM NUMBERS AND FINITE SIZE EFFECTS IN THE GROUNDSTATE 



Let us start with the groundstate properties of the Hamiltonian (1.1) in the strip 



0<a<i. (2.1) 

In the absence of a magnetic field, the groundstate is found in the sector with total spin Sz — and momentum pa 
{po — 0,N — 2n, n even, pQ = n,N = 2n, n odd). We obtain the ground state energies E{Sz,a, A, N) on finite systems 
up to A^ = 30 through a direct Lanczos diagonalization, making use of the translational invariance of Hamiltonian 



(1.1). This reduces the dimension of the Hilbert space approximately by a factor of N. We choose a set of base states 
as proposed in Ref. E^ by Takahashi. This choice results in a real Hamiltonian matrix even for momenta p 7^ 0, tt. 

Hence we obtain the proper groundstate of the model even if its momentum is not p = ot p = n. 

Along the line a = i, A > — i the groundstate is twofold degenerate; the_two states are just the dimer states (pT^), 
(1.3). The proof of this statement follows the arguments of van den BroeUlil. The Hamiltonian 

i/(a = i,A,B = 0) = ^ff(*,i + l,z + 2,A) (2.2) 

i 

is expressed in terms of three spin Hamiltonians: 

Hi^, z + 1, z + 2, A) = \{S+Sr^, + SrS+,, + S+,,Sr^, + Sr^.S^, + S+Sr^, + SrS^,) 

The dimer states ( [l.2[ ), ( |l.3D turn out to be eigenstates of H{i,i + 1,« + 2, A) with eigenvalue eo(A) = — (^ + 4-). 
One can easily prove that this is the lowest eigenvalue of the three spin Hamiltonian for A > — i. However, this is 



not the case for A < — i, where the lowest eigenvalue of the three spin Hamiltonian is given by ei(A) — —^- The 



corresponding eigenstate has all three spins pointing in the same direction (cf. (|l.6|)). 

Let us next turn to the finite size effects of the groundstate energy. Within the strip 0<a<^;A>— i, they turn 
out to be monotonically increasing with N. Finite size effects vanish on the 'dimerline' a — ^i ^ ^ ~2' where the 
groundstate is completely dimerized and degenerate. Right to the dimerline, i.e. for a > i, A > — i, the groundstate 
properties - with respect to its momentum quantum numbers - change and the monotonic behaviour of the finite size 
effects is lost. 

In the presence of a uniform magnetic field B with magnetization M{B) = Sz/N the ground state of the isotropic 
model (A = 1) is found iii^he sector with total spin Sz ~ S. A rule for the momenta ps of these states can be deduced 
from Marshall's sign ruleca: 

Ps^O for 2S + N ^ An 
Ps^TT for 2S + N = 4:n + 2. (2.4) 

This rule has been provenQU to be correct in the unfrustrated case, however it turned out to be valid in a larger 
M-d ependent domain in the a — A— plane. E.g. in the isotropic case (A = 1) we foundEl that the momentum rule 



(2.4) is satisfied for a < ao{M) i.e. below some curve ao(M), which starts at the Majumdar Gosh point 



and ends at 



ao(M = 0,A = l) = i (2.5) 

«o(Af = i,A = l) = i. (2.6) 

The groundstate is degenerate along the curve ao(M, A). The two states differ in their momenta; the first one follows 



(2^). A rule for the momentum of the second state has not yet been found. 

In Fig. 1 we have plotted numerical results of the curves ao{M,A) for various values of the anisotropy parameter 
A = 1.0,0.4,0.1,-0.2. The data for ao{M,A) mark those points in the a — M-plane, where the groundstates with 
energy e{M,ao{M,A),A,N), N = 12, ...,18 are twofold degenerate. Finite size effects of ao{M,A) are visible at 
M = 1/4 and M = 1/6, 1/3 where systems of size N = 12, 16 and N = 12, 18, respectively, are realized. In spite of 
the finite size effects we think that the finite system results shown in Fig. 1 reproduce the qualitative features of the 
curves ao(M, A) in the thcrmodynamical limit: 



• all curves start and end at the points (|2.5|) and (p.6|). 



• aolM, A = 1) has a pronounced maximum around M = 0.2 with rather large finite size effects. For decreasing 
values of A the height of the maximum is reduced and its position is shifted to smaller values of M . 



Beyond the curve ao{M, A) — i.e. for a > ao{M, A) — the groundstate momenta deviate from the rule (2^) and we 
therefore expect a change in the groundstate properties. 

III. THE PHASE DIAGRAM IN THE PRESENCE OF A MAGNETIC FIELD 

In this section we will present numerical results for the groundstate energy per site e{M,a, A, N) — 
jfE{Sz, a, A, TV). We are in particular interested in the changes of the M-dependence of these energies with a and A 
since they indicate a change in the groundstate ordering. The following situations have been found: 

• Ferromagnetic phase: 
Here the free energy 

f{M,a,A)^e{M,a,A)-BM (3.1) 

is minimized by the states \F± >, where all spins are up (+) or down (— ), respectively. It turns out that the 
boundary of the ferromagnetic phase A < A/ (a) is characterized by the degeneracy 

e{M = 0, a, Af{a), N) ^ e{M = i a, A/(a), N) (3.2) 

of the lowest energy eigenvalues in the sectors with Sz ~ and Sz ~ N/2. 



Antiferromagnetic phase: 

The minimum of the free energy is found for Q < M < ^ at 



"'^ >0 (3.4) 



dA/2 
This means that e{M, a, A, N) is monotonically increasing and convex for < M < i The saturating field 



de . „,., 1 

dM 



M=i^BiM^-,a,A), (3.5) 



which is needed to ahgn all spins in the system, can be computed from the one magnon states: 

\p,S. = ^-l>=^Y.e^P-\x>, (3.6) 



2 VN 

where \x > denotes the state with one spin down at site x and all other spins up. The energy of this state is 

N N 

E{S^ = — -l,a,A,N) = cosp + a cos 2p + A(l -«)(—- 1) (3.7) 

and the groundstate energy is found by minimization with respect to p. For < a < 1/4 the minimum is found 
But p — TT and the saturating field is 

B{M=-,a,A) = A{l + a) + {l-a). (3.8) 



For T < a < i the minimum energy (3.7) is found for 



cosp = (3.9) 

4a 



which yields for the saturating field 



B{M ^ -, a, A) ^ Ail + a) + a +—. (3.10) 

2 sa 

The boundary Aa{a) of the antiferromagnetic phase A > Aa{a) is characterized by the condition 

^(M,a,Aa{a),N)\,J^^=Q (3.11) 

i.e. the convexity condition is lost for A < Aa{a). 

• Metamagnetic phase: 
Between the ferromagnetic and the antiferromagnetic phase 

A/(a) < A < Aa{a) (3.12) 

we find a metamagnetic phase, which is characterized be a zero in the second derivative: 

d'^e 
> 0, < M < Mcia, A) (3.13) 

'^ ' (Af,a,A,7V)|M=M. =0 (3.14) 



dAP 

d\ 
MP 

The minimum of the free energy is found for 



-^^B for < M < M,{a, A) (3.15) 

and at 

M=l for B>B, = -^\m=m^. (3.16) 

Therefore, in this metamagnetic phase we have a discontinuity at Bc{a, A) where the magnetization curve jumps 
from M = Mc{a,A) to M — i. Bc{a,A) decreases, if one crosses the metamagnetic phase coming from the 
antiferromagnetic phase and moving towards the ferromagnetic phase. An example will be given below. 
For small magnetic fields < B < B^. the system looks antiferromagnetic, for B > B,, ferromagnetic. 

For the determination of the phase boundaries Af{a) and Aa{a) of the ferromagnetic and antiferromagnetic phase, 
we have first computed the lowest energy densities e(M, a, A, N) as they depend on the magnetization M and the 
parameters a and A. As an example we show in Fig. 2a) the evolution of the M-dependence for A = —0.6, a = 
0.0, ..., 0.40 on a system with A^ = 16 sites. One clearly observes the three phases discussed above. For small values of 
a, e(M, a. A) is monotonically increasing and convex. Here, we are in the antiferromagnetic phase. At a = 0.15 the 
second derivativ e (|3.14| ) vanishes first at Mc = 5. The metamagnetic phase extends from a = 0.15 to a = 0.305, where 
the degeneracy ( |3.2| ) shows up. The behaviour of the corresponding magnetization ciujves can be seen in Fig. 2b). 
We obtain these magnetization curves by applying the method of Bonner and Fishertj to our finite system results. 
For a > 0.305 we then enter the ferrom agn etic phase. The resulting phase diagram in the a — A— plane is shown 



in Fig. 3 . The numerical evaluation of (3^) on finite systems does not show a significant finite size dependence. In 



other words, the determination of the phase boundary Af{a) is well u nder control. 



The determination of the second phase boundary Aa{a) from (3.11) turned out to be much more difficult. We 
numerically calculated e{M = 1/2 — 2/A^, a, A, N) - i.e. the lowest eigenvalue in the sector with two spins flipped - 
on rather large systems with A^ = 20, 30, 40, 50 and looked for a zero in the second derivative: 

eiM = 1 - 1, a, A, TV) + e{M = 1, «, A, A^) - 2e(A/ = i - 1 a. A, A^) = (3.17) 

The resulting Aa(a) suffers under finite size effects particularly in the vicinity of the point a = 0.5, A = —0.5. The 
curve Aa{a) plotted in the phase diagram (Fig. 3) represents the result of (|3.17) for the largest system size A^ = 50. 



The points a = 0, A = — 1 and a = i, A = — i are special in the sense, that the boundaries Aa{a),Af{a) for the 
antiferromagnetic and ferromagnetic phase meet. Therefore, we have no metamagnetic phase between the ferro- and 
antiferromagnetic phase at these points. This can also be clearly seen in Figs. 4a and 4b, where we have plotted the 
M and A dependence of e{M,a = 0, A) and e{M,a — 1/2, A), respectively. These energies turn out to be convex 
for A > —1, (a — 0) and A > — i, (a — 0.5) and concave for A < —1, {a = 0), A < — i, {a — 0.5), respectively. 
At A = —1, a = and A — — i, a = i, the groundstate energies e(M, a, A, A^) are all degenerate with respect to M. 
This feature will be investigated in the next section. 

IV. GROUNDSTATE DEGENERACY AT THE POINTS WHERE THE PHASE BOUNDARIES MEET 



According to ( |3.2[ ) the phase boundary A/ (a) of the ferromagnetic phase is defined by the degeneracy of two 
eigenstates with total spin Sz = and Sz — ^ ^ respectively. At the points A = — 1, a = and A = — i, a = ^ where 
Af{a) — Aa{a), a much larger degeneracy of the groundstate with respect to all values of 6*2 = 0, ..., -y occurs. We 
are going to show now, that the n-magnonstates S-f^{p)"'\F— > — obtained by n-fold application of the rising operator 

S+ip) = Y,e^P'S+ (4.1) 

on the ferromagnetic state \F~ > — are eigenstates of the Hamiltonian 

H{a,A,B)S+{pr\F- >= EpS+ipriF- > (4.2) 

for 

a) a = 0, A = -l, p = 7r (4.3) 

1 2tt 47r 

b) -00 <a<oo, A = --, ^=-3-' ^^"3"' ^^'^'' 



Here 



N 
Ef = — (A + A-a) 



(4.5) 



is the energy of the ferromagnetic state. For the proof of (4.2) we start from the commutation relations 



[H{a, A, B), S+{p)] = 5] (a, ^ S^ S^^^e^^' + b, J^ ^f ^z+,e'"') 
j=i I I 



where 



and 



(4.6) 



ai = -e'P + A, 5i = -1 + A( 



ip 



a2 = a{~e^'P + A) 62 = a(-l + Ae^'P) 



(4.7) 



where 



[[iJ(a, A), S+ip)],S+ip)] = E ^^- E ^;^^i,e'^'' 



ci = 2e*''(A-cosp) 
C2 = 2ae2*P(A-cos2p) 



(4.8) 



(4.9) 



All further commutators with S^ 
\F- > yields: 



(p) vanish identically. Application of (|4.6|) and (fl^) onto the ferromagnetic state 



[H{a,A),S+{p)]\F~ >=-2(A-cosp + a(A-cos2p))|p> 



where \p > is the one-magnon state (p.6[). Similarly one finds 



where 



[[i7(a,A),5+(p)],5+(p)]|^- > =2e'P(A-cosp)|2p,l > 

+2e^*Pa(A-cos2p)|2p,2 > 



\2p,j>=J2e''''%l+J> 



(4.10) 



(4.11) 



(4.12) 



are two magnon states with two spins up at sites I and I + j. Note, that the right hand sides of (4.10) and (4.11) 
vanish for the two cases listed in (p^) and (4.4) , respectively. The first point ( [4. 3D marks the transition from 
antiferromagnetism to fcrromagnctism in the nearest neighbour XX Z model. Indeed, the rising operator >5'+(7r) 
commutes with the Hamiltonian H{a = 0, A = — 1). 

Along the line (4.4), the e ndp oint of the dimcrlinc (a = 0.5, A — — 0.5) is of special interest. Here the eigenvalues ( |l.5| ) 
of the dimerstates (1.2), (L3) are degenerate with the energy (15) of the ferromagnetic states. This also implies that 
the n-magnonstates ( |4.2| ) are groundstates for a = i which explains the degeneracy found in Fig. 4 a fo r a — —0.5. 
In Fig. 5 we have plotted the groundsta te e nergies e(M, a, A = —0.5) on a ring with 18 sites along the line (4^), where 
the degeneracy of the n-magnon states ( |4.2| ) has been proven. For a < 0.5 the groundstate energies e(M, a, A = —0.5) 
are monotonically increasing with M; i.e. the corresponding groundstates cannot be identified with the degenerate 
n-magnonstates (4.2). The groundstate energies e(Af, a, A = —0.5) meet each other for all M at a = 0.5 and stay 
very close together in the interval 0.5 < a < 0.6. This leads to the narrow width of the metamagnetic phase in Fig. 
3 for 0.5 < a < 0.6. For a > 0.6 the quasi-degeneracy with respect to M is lifted again. 



V. DISCUSSION AND CONCLUSION 



Metamagnetism denotes a mixed phase between ferromagnetism and antiferromagnetism, which has been observed 
in various substances like e.g. FexMni-xTiO^. The characteristic signal is a rapid increase (or discontinuity) in 



the magnetization curve, if the external field exceeds a critical value Be- For B > Be the substance is almost fully 
magnetized. In this paper we have shown that the phenomenon of metamagnetism can be observed in the one 
dimensional spin-1/2 XXZ-model with next to nearest neighbour coupling a and anisotropy A. The phase diagram 
in the a — A-plane (Fig. 3) contains three regimes: the antiferroniagnetic one with A > Aa(a), the ferromagnetic one 
with A < A J (a) and the metamagnetic one in between Aa{a) > A > Af(a). The nietamagnctic phase shrinks to zero 
at a = and a = 0.5, where A.a{a) — Af{a). At these points there is a direct transition from antiferromagnetism to 
ferromagnetism and the groundstate turns out to be highly degenerate — namely with respect to Sz — 0, 1, 2, ..., N/2. 
These states can be identified with n-magnon states. 
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Figure Captions 



FIG. 1. The groundstates of the Hamiltonian (1.1) are degenerate along the curves Q;o(Af, A). The numerical data 



points were obtained on finite system calculations with TV = 12, 14, 16, 18 and are shown for A = 1.0, 0.4, 0.1, —0.2. 



FIG. 2. (a) M-dependence of the groundstate energy per site e(M, a. A, N) with A^ = 16, A = —0.6, and a = 
0.0,0.1, ...,0.4. 

(b) Magnetization curves M{B, a, A) for A — —0.6 as they follow from the groundstate energy per site e(M, a, A, N) 
shown in Fig. 2a . In the metamagnetic phase (e.g. at a = 0.2,0.25, A = —0.6) there is a critical field Bc{a,A), 
where the system jumps into the fully magnetic state with M = 1/2. 



FIG. 3. The phase diagramm in the a — A-plane, as it follows from the numerical evaluation of ( |3.2| ) {N — 18) and 



(3.10) (N ~ 50). The metamagnetic domain A/(a) < A < Aa(a) is hatched. 



FIG. 4. Groundstate energies e(M, a, A, iV = 20, 24) along the hues 

(a)a = 0.5, A = 1.0, 0.6,..., -1.0, (b) a = 0.0, A = 1.0, ..., -1.8 . 

The groundstate degeneracy is clearly visible for a — 0.5, A = —0.5 and a = 0, A = —1, where the direct 

transition from antiferromagnetism to fcrromagnctism occurs. 



FIG. 5. a and M dependence of the groundstate energies e(M, A,N — 18) along the hnc A = —0.5. 
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